Universality in entanglement of quasiparticle excitations 
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We show that the entanglement entropy of single quasiparticle excitations of one dimensional 
systems exceeds the ground state entanglement entropy for log(2), if the correlation length of the 
system is finite. For quadratic fermion systems we show that the excess of entanglement is related 
to the number of quasiparticles in the excited state. This observation is confirmed by numerical 
examples of one dimensional quantum many-body systems, including nonintegrable. Tensor network 
methods to describe quasiparticle excitations are discussed. 



PACS numbers: 03.65.Ud, 05.70.Ce,75.10.Jm 

The entanglement is of central interest in quantum in- 
formation theory, quantum computation and plays an 
important role in quantum many-body physics. On one 
hand, it is a valuable resource in quantum computation 
and the reason why quantum algorithms can be more 
powerful than classical ones |T|. On the other hand, it 
makes classical simulation of quantum many-body sys- 
tems more difficult in terms of tensor network methods 
relying on a low degree of entanglement, such as the 
density matrix renormalization group [5] or matrix prod- 
uct states [21 3] where the number of parameters scales 
roughly exponentially with the entanglement entropy [5] . 
The entanglement of quantum many-body systems at 
zero temperature has been extensively studied and, in 
particular for one dimensional lattice systems, it was 
shown that the ground state is much less entangled than 
generic quantum states, which makes a classical simula- 
tion tractable [5] , and on the other hand provides a signa- 
ture of quantum critical behavior [7 where the entangle- 
ment entropy violates the area law [5], connected to the 
central charge of the corresponding quantum field theory 
[SHTTj. Much less, however, is known about the entan- 
glement of excited states. It was shown by Alba, Fagotti 
and Calabrese [12] that entanglement entropies of excited 
states of integrable systems appear in a band-like struc- 
ture and their ratios with the ground state entanglement 
entropy scale either extensively or logarithmically with 
the block size. For a critical Ising model, an exact scaling 
relation was obtained by Berganza, Alcaraz, and Sierra 
[F3] using a field theoretical approach. Still, no indication 
has been given that the banded structure of the excess 
of entanglement, as the difference from the ground state 
entanglement entropy is called in [13| . comes from the 
quasiparticle picture at low temperatures where a quasi 
particle causes an excess of log(2) to the entanglement 
entropy as compared to the ground state. This is what 
we show in the present Letter, assuming a translation 
invariance, a finite correlation length and a validity of a 
quasiparticle description at low energies. For quadratic 
models we also show that in the thermodynamic limit, 
the entanglement entropy of fc-quasiparticle excitation 



exceeds the entanglement entropy of the ground state 
for fclog(2). The observation is supported by numerical 
simulations of spin chains: the Heisenberg XY model, the 
Ising model in a transverse field and the Ising model in 
a mixed field, the last being a nonintegrable model [14] . 
Finally, we show that the excess of entanglement for the 
Bloch-like Ansatz for excited states [15] is exactly equal 
to log(2) and is as such well suited to describe single 
particle excitations. 

Let us consider a one-dimensional quantum lattice sys- 
tem of length n described by a Hamiltonian H which 
is invariant to spatial reflections, [H, R] = 0, and con- 
serves a parity of the system, [H, P] =0. In a sym- 
metric bipartition of the system into two parts of length 
L = n/2, any pure state can be parametrized by a matrix 
ip as \ip) = J2i j V'ijl*) ® |j) where indices i,j run over 
configuration states of the subsystems. The Schmidt de- 
composition reads \tp) = J2k ^(E* Oi.fcl*))® (Ej Pj,k\j)) 
where -0 = cta/3 H . We denote the ground state as ft 
and assume it is nondegenerate and thus reflection sym- 
metric and of an even parity. In this notation, it is 
evident that the reduced density matrix of the system, 
tr r i g ht[|^)(^|], is simply flfl H and the entanglement en- 
tropy, defined as the von Neumann entropy of the re- 
duced density matrix, reads S [n] = -tr[flfl H \og(flfl H )]. 
In the following we shall be interested in the entangle- 
ment of single particle excitations |$) = a"t|0) where X 
is a quasiparticle annihilation operator fulfilling canoni- 
cal anticommutation relations (CAR) {A,A^} = 1 and 
{X,X} = and flipping the parity, {X,P} = 0. We 
shall call AS = S(Q) — 5(f2) the excess of entanglement, 
as introduced in |13) . Like the ground state, the ex- 
cited states must also obey the symmetries of the sys- 
tem, being either symmetric or antisymmetric to reflec- 
tions, R\$) — ±|$), and of an odd parity. Because of 
these symmetry constraints, there exist maps Xi and 
X r , acting only on the left /right part of the system, 
such that |$) = ^{X\ +Xt)|0), and {X Lr ,P} = 0. 

They anticommute {Xi,X r } = {Xi,X}.} = and thus 
{X ltr ,xl r } = |1 and {Xi,^} = -{X r ,X r }. 
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Let us now assume a finite correlation length £, such 
that \{£ixj) - (xi){xj}\ < Ce^^^/t f or any two local 
operators Xi, xj at sites Although the quasiparticles 
are not localized, simultaneous presence in both halves 
of the systems gets unlikely for n — > oo, 



|<n|x,xt|n)| < 



const 



= 0(n" 1 ) (1) 



which results from X = ^= YH^iLfj^j where {c,} 

are real particle annihilation operators at sites {j} and /j 
and gj are bounded. The actual quasipartices might be 
of a different form but they still behave like real particles 
and this is a valid assumption. This bound also results 



in \\X l Q.\\* F = {Q\X}X l \n) 



-(n\xjx r \n) = o(r 



that = tr g |Q)(fi| since |f2)(f2| is separable in 

the space of quasiparticles H = (£) q H q being the 
Fermi sea). Let us now orthogonalize some H q with 
respect to H r and write the result as H q > < Hi such 
that H = Hi' <E> Hq' ® H r - Obviously, pv+q> and py 
are reduced density matrices for the ground state and 
for the quasiparticle excitation, respectively. Using an 
inequality S(pi>) < S(pi>+ q >) + S(p q >) pQ we arrive at 
S [<f q ] < glto] + log ( 2 ) because p q > is of rank 2 and thus 
S{p q .) < log(2). 

Quadratic systems. An example of systems where the 
quasiparticle picture is not only valid but exact, are in- 
tegrable systems, mappable to free fermions, such as a 
quantum Ising chain in a transverse magnetic field. We 
consider a quadratic Hamiltonian in terms of (hcrmi- 



The reduced density matrix p for the excited state $ = tian) Majorana operators {wi,Wj} — 2Sij on n sites, 
J=($i+$ r ) where $ ; = XfQ and $ r = OX*, now reads R = £ 2 " =1 w.H^Wj, where H = H H = -H T . The 
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The first two terms describe states where either of the 
two parts is excited and have the same spectrum be- 
cause of the reflection symmetry. We will show that 
they are distinguishable whereas the mixed terms van- 
ish in the thermodynamic limit. Let us first fix the 
notation, = Ylj <jj(x x), \\x\\\ = Ylj <J j( x ) ancl 

II^Hoo = (Ti(x) where a\(x) > <i2{x) > ... are sin- 
gular values of x. The distinguishability follows from 
tr[$,$f* r *f] < WnWlWX&X*^ where ||fi|| F = 1 

H\ II V.DII2 _ r>f m -l\ 



and \\x,nx: 



' X r X, 



to ||^0||| = Oin- 1 ) from |T]). We have used tr[pv] < 
I |ootr[zy] < tr[/x]tr[i/] for p, v > (see e.g. [IB]) and 
XrXi 1 < I =>■ 1 1 X r I loo < 1- The trace norm of the 



mixed term <J> ; $f = XfQX^Q H 



-X/*JQQQ ff also 



vanishes, \\X^ XiVtVt H \^ < {{{X^ X t ) 2 (nn 



H\2||l/2 



< 



i oo 



\Xtfl\ 



\M\i 



< 



0{rr x l 2 ). We have used 
for p, v > and t > 1 [16] . It is 
easy to show that the trace of the mixed norms vanishes 
as 0(n" 1 ), again due to (nj). Therefore, 



ImVIIi 



p={ l -u?tm H u l )@( l - 



u ?n H nu r ) + p' 



where \ \p'\\i = 0(n -1 / 2 ) and u^ r are isometries for which 
uf r Ur y i = 0. Because p' vanishes and tr[u H QQ H u] < tr[p] 
for any isometry u and p > 0, whereas on the other hand 
tr[p] = 1, we conclude aj(u^VlVl H ui) — aj(u^fl H ^lu r ) = 
<jj(QQ H ) + 0(n -1 / 2 ). Because of the continuity of von 
Neumann entropy, \S(p) — S(p')\ < \\p — p'||i, the ex- 
cess of entanglement for the excited state $ approaches 
log(2) with corrections of the order of 0(n -1 / 2 ). We will 
learn from numerical simulations that this bound is not 
optimal and the corrections are actually of the order of 

The upper bound AS < log(2) is given by splitting 
the total Hilbert space H — Hi <8> H r and realizing 



uij are chosen such that {w™ 1 • ■ • w™£ L ;rrij e {0,1}} 
only act on the Hilbert space of the first L = n/2 sites 
for which they form a complete basis. We diagonalize 
H = VAV H - V*AV T with a 2n x n isometry V for 
which V H V — I and V T V = 0, and define quasiparti- 
cle annihilators = J3i=i ^fc^i) so ^hat we obtain a 
free fermion model H = J2 k EftStSfc. As shown in [HH2], a 
2Lx2L correlation matrix {ip\wiWj |V) = ^,i+irg where 
is a real skew symmetric matrix with eigenvalues 
±iv^ , determines the eigenvalues of the reduced density 



|Xjn||> = 0{n 1 ), due with m 



matrix p^ = tr righ t | i>) (ip \ as p [ m 



l+rrij v 



±1. The entanglement entropy of ip reads 



SW = -tr(^±f^log^±f^). Let us denote the upper 
part of V, referring to sites {1, . . . , L}, as u E C 2Lxn , 
where we also reshuffle the columns u k in a way that 

hShi, because 



J*2fe-1 mi 



of the reflection symmetry. We use a - b = a b. For the 
ground state, it is easy to check that (wiWj) = 2[uu 



Hi 



J'-./ 



and rl n l = i(u*u T - uu H ). Bases {\/2u 2k _ 1 } and 
{\/2w2fe} are orthonormal which allows us to write 
as a sum of L rank- 2 terms 



r [fi] 



L 

E 

fe=i 



X 2 fc-i for Xk = 2i(v* k ® vl - v k ® ). 



For an excited state |$ K ) = &^|^), it is easy to see that 
r [*»] = r [o] _ Xk = X2fc _ i5 assuming /t = 2/t' 1, 

which is obviously rank deficient with dimkerr^l = 2. 
This makes the spectrum of the reduced density matrix 
double degenerate. The same reasoning is applied to 
^L'-i^'-il Q ) ( or f>L*>2 M 'l n »> resulting in a four times 
degenerate spectrum of the reduced density operator. Let 
us now show that AS approaches log(2). Using the Tay- 
lor expansion, we write 

oo 

AS = 5>j(2j - l)}-Hr[(irM)V - (tTM - IXk ) 2 >]2- 1 

3 = 1 
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FIG. 1. Excess of entanglement AS for the tilted Ising chain 
on n = 100, 200, 500 (top to bottom; green, blue, red) sites for 
reflection symmetric (plus) and antisymmetric (cross) states. 
We have connected points where |A5/log(2) — 1| < 0.1. The 
inset shows AS vs energy for n = 200. 



where tr[(ixfe) 2j ] < 2. If (Vk,l)u k ■ u ; =0, then the 
trace in above expression equals 2, since XkXi = dkiXk? 
and thus AS = log(2). Using b k . L = ^ J2i=i u i,k^i 
we see that this is indeed the case for n — > oo where 

Ei\u k -m? = (n\bl L bk-x\n) = o^- 1 ) due to Q. 

Hence, the excess of entanglement for a single quasipara- 
cle excitation equals AS = log(2) and for fc-quasiparticle 
excitation it equals fclog(2) (if all excited quasiparticlcs 
belong to the same reflection symmetry class). 

Results. From a very generic setting we have seen 
that the entanglement entropy of excited states is con- 
nected with the number of quasiparticles which we shall 
now confirm by numerical simulations. First we con- 
sider an antiferromagnetic Ising chain in a mixed trans- 
verse and longitudinal magnetic field, the so-called tilted 
Ising chain, described by a Hamiltonian operator H = 
J2j ( cr j' T j+i + a j + a j ) where we have chosen a non- 
critical (gapped) regime featuring strong non-integrable 
effects. The model is not exactly solvable and we rely 
on the density matrix renormalization group [2] to calcu- 
late the first few lowest energy excited states. To check 
the validity of the results, we increase the bond dimen- 
sion and optimize the results by a variational optimiza- 
tion technique |17) . Eventually we settle on D = 100. 
The states are rotated in a way that they are eigen- 
states of the reflection operator. We calculate the en- 
tanglement entropies for M = 50 lowest energy excited 
states and various chain lengths, n — 50, 100, 200, 500. 
From Fig. [T] we observe that the excess of entanglement 
approaches AS — > log(2) for the lowest exited states 
which is consistent with single quasiparticle excitations; 
reflection symmetric and antisymmetric excitations dis- 
play distinct behavior. We observe another clustering at 
AS ~ 2 log(2) which corresponds to two-quasiparticle ex- 
citations. Note that for n = 500, all states that we calcu- 



FIG. 2. Excess of entanglement AS* of single quasiparticle 
excitations fejjfi) of the XY(7 = 0.5, h = 0.9) model for n = 
128,256,512 (top to bottom; green, blue, red) and reflection 
symmetric (+, plus) and antisymmetric (— , cross) states. The 
inset shows the scaling of AS with n at phase n/2 for .R* 
(plus/cross). 



late using the DMRG, are single particle excitations since 
the two-particle excitations only start at energy approx- 
imately 0.7 (i.e. two lowest single particle excitations). 
Above this threshold, we also observe reflection invariant 
states with irregular behavior which, we conclude, must 
be bound states. 

We do a systematic study of excitations for an exactly 
solvable Heisenberg XY model, described by a Hamilto- 
nian H = J2j ^f L<J j a j+i + ^f L tftf+i+ HZj ha j where 
we choose a particular point in the parameter space 
(7, h) — (0.5,0.9); any point with /i / 1 and 7 / 
gives similar results. In Fig. [2] we show the excess of 
entanglement of all single particle excitations b\\£l) for 
chain lengths n = 128, 256, 512 and again observe that it 
approaches AS —> log(2). The scaling analysis confirms 
that log(2) — AS scales as n _1 as shown in the inset. 

Three-particle excitations are studied on the case of 
an Ising chain in a transverse field (i.e. the XY model 
with (7, h) = (1,2)). We choose excitations of a form 
|$) = blty, bf, |f2) where k\ and ki are fixed and corre- 
spond to reflection antisymmetric excitations at quasi- 
momenta tt/4 and 37r/4. We observe in Fig. [3] that 
AS* — > 3 log(2) for n ~ > 00 except for two peaks at i ~ k% 
and i ~ &2 where two quasiparticles can form a bound 
state and the excited state looks like a two-particle exci- 
tation with the excess of roughly, but more than 2 log(2), 
as the numerical results suggest. 

Discussion. The main message of this Letter, that the 
entanglement entropy of single particle excitations ex- 
ceeds the the ground state for log(2), can be used to jus- 
tify the Ansatz for single particle excitations |15j. They 
are based on a homogeneous matrix product state for 
the ground state described by a tensor A s e £DxD ag 
l^> = J2s v i(--- AS - lAS ° ASl ---) v r\s), see [IB], and are 
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FIG. 3. Excess of entanglement AS* of bjb^^^yjfi) for the 
Ising model (h — 2) with n — 128, 256, 512 (top to bottom; 
green, blue, red) for reflection symmetric (R + , plus) and an- 
tisymmetric (R~ , cross) states. 



given as a sum of matrix product states where all ten- 
sors are equal to A a except for a site j where some other 
tensor B s is used, together with a phase factor e KJ for 
momentum n. Formally, the Ansatz reads |<& K (-B)) = 
l z r 1/2 'Emez-s e lKm f m v\{- ■ ■ A s -^B Sa A s ^ ■ ■ ■ )v r \s). For 
each k, there exist D 2 (d — 1) (d = 2 for spin-1/2) or- 
thonormal states !$«(£)}, such that (Q,\$ K (B)) = 0. The 
relevant quantities in infinite matrix product states are 
the left and the right eigenvector of a transfer matrix 
E = ^ S A S ® A , corresponding to the largest eigen- 
value, which can be reshaped to hermitian positive def- 
inite D x D matrices l,r. An overlap of semi- infinite 
states with respect to a bipartite splitting, fl = a(3 T 
in the notation of the first part of the Letter, is given 
exactly by these matrices as a H a — I and f3 H (3 = r T , 
and so is the reduced density matrix which reads p = 
a[3 T f3*a H — 6iE.a H where 3 = L H rL and a = aL~ H 
such that a H a — I. We have used a Cholesky decom- 
position I = LL H . The excitations can be written as 
$ = -^{4>P T + oi9 T ) where the first/second part con- 
tains all excitations on the left /right part of the chain, 
ft is easy to show that also cj) H <fi — I and 6 H 9 = r T 
whereas 9 H a = <fi H (3 = (and <f> H a = by construction) 
which means that 4> is orthogonalized by the same trans- 
formation as a, that is cf> = 4>L~ H . The reduced density 
operator reads p — ^{<f)r(f) H + ara H ) and we immedi- 
ately obtain p = l;((f)'E(f) H + aE,a H ) where and a are 
isometries with orthogonal support, a H <f) = 0, and thus 
p ~ !(3©H). The entanglement entropy of the Ansatz is 
equal to Si*! = S™ + log(2), irrespectively of the bond 
dimension D, which we have shown to agree with single 
quasiparticle excitations in the thermodynamic limit. By 



this we provide additional support that the Ansatz in |15j 
is a good Ansatz for describing single particle excitations. 
Unfortunately, however, the entanglement entropy of al- 
most all other excited states (except e.g. some bound 
states) goes beyond this limit; in this case a generalized 
Ansatz should be used where not one but two or more 
tensors in A 3 the ground state are replaced by some B s . 

Conclusion. We have studied the entanglement en- 
tropy of single quasiparticle excitations and shown that 
it is related to the ground state entanglement entropy 
with a contribution of log(2) for each quasiparticle in 
the thermodynamic limit, assuming a finite correlation 
length and a valid quasiparticle description at low ener- 
gies. We have provided numerical evidence for this uni- 
versal behavior which suggests a possibility to classify 
excited states on the basis of the entanglement. 
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